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Abstract: Sparticle mass hierarchies will play an important role in the type of signatures 
that will be visible at the Large Hadron Collider. We analyze these hierarchies for the four 
lightest sparticles for a general class of supergravity unified models including nonuniversalities 
in the soft breaking sector. It is shown that out of nearly 10 4 possibilities of sparticle mass 
hierarchies, only a small number survives the rigorous constraints of radiative electroweak 
symmetry breaking, relic density and other experimental constraints. The signature space 
of these mass patterns at the Large Hadron Collider is investigated using a large set of 
final states including multi-leptonic states, hadronically decaying rs, tagged b jets and other 
hadronic jets. In all, we analyze more than 40 such lepton plus jet and missing energy 
signatures along with several kinematical signatures such as missing transverse momentum, 
effective mass, and invariant mass distributions of final state observables. It is shown that 
a composite analysis can produce significant discrimination among sparticle mass patterns 
allowing for a possible identification of the source of soft breaking. While the analysis given 
is for supergravity models, the techniques based on mass pattern analysis are applicable to 
wide class of models including string and brane models. 
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1. Introduction 

Supersymmtery (SUSY) remains a leading candidate to describe new physics beyond that 
of the Standard Model (SM). Recently, an approach for identifying supersymmetric particles 
(sparticles) was proposed involving sparticle mass hierarchies, or sparticle mass patterns. 
Such patterns could yield distinct identifiable signatures at the Fermilab's Tevatron and at the 
CERN Large Hadron Collider (LHC)(l|, ||. At the same time, the hierarchical mass patterns 
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are model dependent and the determinations of such patterns could be helpful in extrapolating 
the data back to the theoretical model. This new approach has been investigated within the 
framework of gravity mediated breaking of supersymmetry^, |], |5| and specifically within the 
minimal supergravity grand unified model, the mSUGRA model || (for a review see Q) with 
sparticle mass ranges that lie within reach of the present colliders (for a review of recent search 
strategies see 0). The analysis of [|| was a rather brief introduction to the technique. Here 
we carry out a more in depth analysis within models with both universal and nonuniversal 



soft supersymmetry breaking J|, ||, 1C[. Thus, in the minimal supersymmetric extension of 
the Standard Model (MSSM) there are 32 supersymmetric particles. We list them here to set 
notation. There are 4 Higgs boson states, of which three (h, H, A) are neutral, the first two 
being CP even and the third CP odd, and one charged Higgs H^. In the gaugino-Higgsino 
sector there are two charged mass eigenstates (charginos) Xi=i2' f° ur charge neutral states 
(neutralinos) Xi=i 41 an d the gluino g. In the sfermion sector, before diagonalization, there are 
9 scalar leptons (sleptons) which are superpartners of the leptons with left and right chirality 
and are denoted as: {eL,R, fiL,R,TL t R,D eL ,D^ L ,D TL }. Finally there are 12 squarks which are 
the superpartners of the quarks and are represented by: {ul,r, cl,r, ^l,r, d^^, sz^r, 6l,,r}. 
Mass diagonal slepton and squark states will in general be mixtures of L, R states. 

If the 32 masses are treated as essentially all independent, aside from sum rules (for a 
pedagogical analysis on sum rules in the context of unification and RG analysis see |12|) 
on the Higgs, sfermions, chargino and neutralino masses, then without imposition of any 
phenomenological constraints, the number of hierarchical patterns for the sparticles could be 
as many as O(10 28 ) or larger. This represents a mini landscape in a loose way reminiscent of 
the string landscape (which, however, is much larger with as many as O(10 1000 ) possibilities). 
[Here we refer to the landscape of mass hierarchies and not to the landscape of vacua as is 
the case when one talks of a string landscape. For the string case the landscape consists of 
a countably discrete set, while for the case considered here, since the parameters can vary 
continuously, the landscape of vacua is indeed much larger. However, our focus will be the 
landscape of mass hierarchies.] Now, the number of possibilities can be reduced by very 
significant amounts in supergravity models with the imposition of the constraints of radiative 
electroweak symmetry breaking (REWSB) 1 , and other phenomenological constraints. This 
was precisely what was accomplished in the analysis of Q . The analysis of Ref. |l], focused 
on the mass hierarchies for the first four lightest sparticles, and found the residual number 
of hierarchies to be 22 in mSUGRA. Here, the possible signatures from some of the patterns 
were also discussed along with the prospects for direct detection of dark matter within various 
mass hierarchies. 

The phenomenology of supergravity (SUGRA) models has been discussed since their 
inception and there exists a considerable amount of literature regarding the implications of 



SUGRA (for early works see [11, 13], for more recent works see 14 , 1£ , |lq| , for works with 



1 EWSB can be realized non-radiatively for certain choices of parameters in the presence of nonuniversal- 
ities in the Higgs sector. Since in this analysis boundary conditions have been imposed at the GUT scale 
and RGEs have been used to obtain the low energy physics, we will retain this terminology in the subsequent 
descriptions of EWSB. 
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nonuniversalities see 17], and for works with hierarchical breaking and with U(l) gauge 
extensions see [18, 2C]). While many analyses of the mSUGRA parameter space have 
been limited to the case of vanishing trilinear couplings, several recent works 22, ^3|, 24. 
O, 16, 25, 26] have appeared relaxing this assumption, and new portions of the parameter 



space have been found consistent with all known experimental constraints on the model. 



In this paper we give a more exhaustive analysis of sparticle mass hierarchies for SUGRA 
models including nonuniversalities and also carry out a more detailed analysis of the signatures 
arising from these patterns. We further focus on ways in which patterns can be discriminated 
from each other using the relevant distinctive features of the signature space. It is found 
that for some model points one encounters the phenomenon where two distinct points in the 
parameter space of soft breaking may yield the same signatures within a 2a error bar. We also 
discuss in this paper how such signature degeneracies can sometimes be lifted by an increased 
integrated luminosity. Finally, we discuss the issue of how well the soft parameters mo and 
m i/2 (where m-o, m^j are the mass parameters in mSUGRA models defined in Sec. (2)) may 
eventually be determined at the LHC which allows one to obtain an estimate on the resolution 
of these parameters using optimal LHC luminosities. 



The outline of the rest of the paper is as follows: In Sec.(||) we give a discussion of the 
sparticle landscape. Specifically, the landscape for the 4 lightest sparticles (in addition to the 
lightest Higgs boson) for the mSUGRA case is discussed in Sec. (2.1) and the landscape for the 
4 lightest sparticles for the nonuniversal SUGRA case is discussed in Sec. (]2.2|) . This includes 
cases with nonuniversalities in the Higgs sector, nonuniversalities in the third generation sec- 
tor, and nonuniversalities in the gaugino sector. In Sec. (|2.3| ), we also discuss the possible 
number of mass patterns that can arise for the above case, as well as for the case when all 
32 sparticle masses are taken into account. An analysis of the patterns and their origin in 
the space of soft breaking parameters is given in Sec.(|3|). An analysis of the benchmarks for 
the landscape of 4 sparticle patterns is given in Sec.( |3.2|) . Sec.(||) is devoted to the discus- 
sion of the sparticle signatures at the LHC. In Sec. ( |4.1| ) we give a discussion of the various 
SUSY tools that are utilized in this analysis. We discuss technical details of the analysis of 
the LHC signatures we have investigated in Sec.(|4,2j). We then move on to discuss how one 
can distinguish sparticle mass patterns arising in mSUGRA in Sec,( |4.3j ), and sparticle mass 
patterns in SUGRA with nonuniversalities in Sec. (|4.4| ). The trileptonic signal as a tool to 
distinguish patterns is discussed in Sec. fl4,5|) . We utilize both event counting signatures and 
kinematical signatures in our analysis, the latter being discussed in Sec. (|4.6|) . In Sec.( |4.7D , a 
method for distinguishing patterns utilizing a large set of signatures is also given. We discuss 
the signature space degeneracy among different models and how to lift it in Sec.( |5.lD , and 
then we generalize our analysis to investigate the resolving power of the LHC with regards 
to its ability to probe the soft parameter space in Sec. (5. 2). Conclusions are given in Sec.(|6|). 
Some of our longer tables have been relegated to the Appendix. 
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2. The Sparticle Landscape 



The analysis proceeds by specifying the model input parameters at the GUT scale, Mq ~ 
2 x 10 16 GeV, (no flavor mixing is allowed at the GUT scale) and using the renormalization 
group equations (RGEs) to predict the sparticle masses and mixing angles at the electroweak 



scale. The RGE code used to obtain the mass spectrum is SuSpect 2.34 [27], which is the 
default RGE calculator in MicrOMEGAs version 2.0.7 [^]. We have also investigated other 
RGE programs including ISASUGRA/ISAJET @, SPheno Q and SOFTSUSY @. We 



have cross checked our analysis using different codes and find no significant disagreement in 
most regions of the parameter space. The largest sensitivity appears to arise for the case of 
large tan f3 and the analysis is also quite sensitive to the running bottom mass and to the top 
pole mass (we take m^ s (m;,) = 4.23 GeV and mt(pole) = 170.9 GeV in this analysis). Such 
sensitivities and their implications for the analysis of relic density calculations are well known 
in the literature and a detailed comparison for various codes can be found in Refs. (]33[|, 

a, n, ©). 

Below we give the relevant constraints from collider and astrophysical data that are 
applied throughout the analysis unless stated otherwise. 

1. WMAP 3 year data: The lightest R-Parity odd super symmetric particle (LSP) is as- 
sumed charge neutral. The constraint on the relic abundance of dark matter under the 
assumption the relic abundance of neutralinos is the dominant component places the 
bound: 0.0855 < n~oh 2 < 0.1189 (2a) |7j- 

2. As is well known sparticle loop exchanges make a contribution to the FCNC process 
b — ► 57 which is of the same order as the Standard Model contributions (for an update 
of SUSY contributions see [^]). The experimental limits on b — > 57 impose severe 
constraints on the SUSY parameter space and we use here the constraints from the 
Heavy Flavor Averaging Group (HFAG) fl| along with the BABAR, Belle and CLEO 
experimental results: Br(b — ► 57) = (355 =t 24^ ± 3) x 10~ 6 . A new estimate of 
Br(B Y s7 ) at 0(a 2 s ) gives Q Br(b -> s 7 ) = (3.15 ± 0.23) x 10~ 4 which moves 
the previous SM mean value of 3.6 x 10~ 4 a bit lower. In order to accommodate this 
recent analysis on the SM mean, as well as the previous analysis, we have taken a 
wider 3.5<r error corridor around the HFAG value in our numerical analysis. The total 
Br(B —* X s7 ) including the sum of SM and SUSY contributions are constrained by this 
corridor. With a 2a corridor, while some of the allowed points in our analysis will be 
eliminated, the main results of our pattern analysis remain unchanged. 

3. The process B s —* ^ + /U~ can become significant for large tan/3 since the decay has 

a leading tan 6 [3 [[ll]] dependence and thus large tan j3 could be constrained by the 

experimental limit Br(B s -► n + fi~) < 1.5 x 10" 7 (90% CL), 2.0 x 10~ 7 (95% CL) |42]. 

-7 



This limit has just recently been updated [43] and gives Br(B s — > fi + fi ) < 1.2 x 10 



(95% CL). Preliminary analyses [44] have reported the possibility of even more stringent 



constraints by a factor of 10. We take a more conservative approach in this analysis 
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and allow model points subject to the bound Br(B s — ► fi + fi ) < 9 x 10 6 (for a review 
see @U). 

Additionally, we also impose a lower limit on the lightest CP even Higgs boson mass. 
For the Standard Model like Higgs boson this limit is ~ 114.4 GeV [fig], while a limit of 
108.2 GeV at 95% CL is set on the production of an invisibly decaying Standard Model 



like Higgs by OPAL |47|]. For the MSSM we take the constraint to be m h > 100 GeV. A 
relaxation of the light Higgs mass constraint by 8 - 10 GeV affects mainly the analysis 
of SUGRA models where the stop mass can be light. However, light stops are possible 
even with the strictest imposition of the LEP bounds on the SM Higgs Boson. We 
take the other sparticle mass constraints to be m~± > 104.5 GeV jp| for the lighter 
chargino, ra^ > 101.5 GeV for the lighter stop, and m^ > 98.8 GeV for the lighter 
stau. 

In addition to the above one may also consider the constraints from the anomalous magnetic 
moment of the muon. It is known that the super symmetric electroweak corrections to — 2 
can be as large or larger than the Standard Model electroweak corrections^]. The implica- 
tions of recent experimental data has been discussed in several works (see, e.g.]50[|). As in 
p5}|, here we use a rather conservative bound —11.4 x 10~ 10 < g^, — 2 < 9.4 x 10~ 9 . 

2.1 The mSUGRA landscape for the 4 lightest sparticles 

One mSUGRA model is a point in a 4 dimensional parameter space spanned by mo, m^rt, 
Aq, tan/3, and the sign of /u, where tuq is the universal scalar mass, is the universal 

gaugino mass, Aq is the universal trilinear coupling, tan (3 is the ratio of the two Higgs VEVs 
in the MSSM, and /i is the Higgs mixing parameter that enters via the term 11H1H2 in 
the superpotential. Typically scans of the parameter space are done by taking a vanishing 
trilinear coupling, and/or by looking at fixed values of tan (3 while varying (mo, m^). In this 
work we carry out a random scan in the 4-D input parameter space for fixed signs of fi with 
Monte Carlo simulations using flat priors under the following ranges of the input parameters 

0<m <4TeV, < m 1/2 < 2 TeV \A /m \ < 10, 1< tan/3 < 60. (2.1) 

Since SUGRA models with \x > are favored by the experimental constraints much of the 
analysis presented here focuses on this case. Specifically for the \i > mSUGRA case, we 
perform a scan of the parameter space with a total of 2 x 10 6 trial parameter points. We 
delineate the patterns that emerge for the first four lightest sparticles. Here we find that at 
least sixteen hierarchical mass patterns emerge which are labeled as mSPs (minimal SUGRA 
Pattern). These mSPs can be generally classified according to the type of particle which is 
next heavier than the LSP, and we find four classes of patterns in mSUGRA: the chargino 
patterns (CP), the stau patterns (SUP), the stop patterns (SOP), and the Higgs patterns 
(HP), as exhibited below 

1. Chargino patterns (CP) : mSPl, mSP2, mSP3, mSP4 
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Table 1: Hierarchical mass patterns for the four lightest sparticles in mSUGRA when /i < and 
H > 0. The patterns can be classified according to the next to the lightest sparticle. For the mSUGRA 
analysis the next to the lightest sparticle is found to be either a chargino, a stau, a stop, a CP even/odd 
Higgs, or the next lightest neutralino y°- The notation A/H stands for either A or H. In mSP14- 
mSP16 it is possible that the Higgses become lighter than the LSP. Y stands for appearance of the 
pattern for the sub case. 



2. Stau patterns (SUP) : mSP5, mSP6, mSP7, mSP8, mSP9, mSPIO 

3. Stop patterns (SOP) : mSPll, mSP12, mSP13 

4. Higgs patterns (HP) : mSP14, mSP15, mSP16. 

The hierarchical mass patterns mSPl-mSP16 are defined in Table (|l]). We note that the 



pattern mSP7 appears in the analyses of [51, 52, 53|. We also performed a similar scan for 



the mSUGRA with [i < case using the Monte Carlo simulation with flat priors and the 
same parameter ranges as specified in Eq.([2.1|). Most of the mSP patterns that appear in the 
fi > case also appear in the [i < case (see Table (|l])). However, in addition one finds new 
patterns shown below 
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mSUGRA^>0 
mSUGRA ii<0 



Landscape of Sparticle 
Mass Hierarchies 




10 11 12 13 14 15 16 17 18 18 20 21 22 



mSUGRA Pattern (mSP) 

Figure 1: Distribution of the surviving hierarchical mass patterns in the landscape for the mSUGRA 
model with fi > (light) and jj. < (dark), under various constraints as discussed in the text. 



1. Stau patterns (SUP) : mSP17, mSP18, mSP19 

2. Stop patterns (SOP) : mSP20, mSP21 

3. Neutralino patterns (NP) : mSP22. 

We note that the analysis of Ref.|34|] has a sparticle spectrum which corresponds to mSPll 
and contains light stops. Light stops have also been discussed recently in |55|, p5|| . 

While the earlier works which advocated benchmark points and slopes made good progress 
in systematizing the search for supersymmetry, we find that they do not cover the more broad 
set of possible mass hierarchies we discuss here. That is, many of the mSP patterns do not 
appear in the earlier works that advocated benchmark points for SUSY searches. For example 
the Snowmass mSUGRA points (labeled SPS) [^6| and the Post-WMAP benchmark points 



of [57], make up only a small fraction of the possible mass hierarchies listed in Table (||). 
The CMS benchmarks classified as Low Mass (LM) and High Mass (HM) p^] (for a recent 
review see ^0[) does a good job covering the mSPl pattern which appears as the most 
dominant pattern in our analysis, but there are no Higgs patterns or stop patterns discussed 
in the CMS benchmarks as well as in SPS or in Post-WMAP benchmarks. We exhibit the 
mapping of mSPs with other benchmarks points in a tabular form in Table (^). 

In Fig. (|l|) we give the relative distribution of these hierarchies found in our Monte Carlo 
scan. The most common patterns found are CPs and SUPs, especially mSPl and mSP5. 
However there exists a significant region of the parameter space where SOPs and HPs can be 
realized. The percentages of occurrence of the various patterns in the mSUGRA landscape for 
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Table 2: Mapping between the mSPs and the Snowmass, Post-WMAP3, and CMS benchmark points. 
The points B 1 = LM1, I' = LM2, C = LM6. HM1 in SuSpect has m^o > m ?1 , but this is not the 
case for ISAJET, SPheno, and SOFTSUSY. Among the CMS benchmarks, only LM1, LM2, LM6, and 
HM1, HM2 are capable of giving the correct relic density. Thus the mapping above applies only to the 
mass pattern, while all of our mSP and NUSP benchmark points satisfy the relic density constraints 
from MicrOMEGAs with SuSpect. The CMS test points do a better job of representing mSPl which 
is the dominant pattern found in our analysis. There are no HP test points or SOP test points in any 
of the previous works. 



both fi positive and [i negative are exhibited in Fig. ([!]). The analysis of Fig. (|l]) shows that 
the chargino patterns (CP) are the most dominant patterns, followed by the stau patterns 
(SUP), the stop patterns (SOP), and the Higgs patterns (HP). In contrast, most emphasis 
in the literature, specifically in the context of relic density analysis, has focused on the stau 
patterns, with much less attention on other patterns. Specifically the Higgs patterns have 
hardly been investigated or discussed. The exceptions to this, in the context of the Higgs 
patterns, are the more recent works of Refs. Q |2|, and similar mass ranges for the Higgs 
bosons have been studied in |6l|] (see also fl6^1 )- 

2.2 The landscape of the 4 lightest sparticles in NUSUGRA 

Next we discuss the landscape of the 4 lightest sparticles for the case of nonuniversal super- 
gravity models. Here we consider nonuniversalities in the Higgs sector (NUH), in the third 
generation sector (NU3), and in the gaugino sector (NUG). Such nonuniversalities appear 
quite naturally in supergravity models with a non-minimal Kahler potential, and in string 
and D-Brane models. The parametrization of the nonuniversalities is given by 

NUH : M Hu = mo(l + S Hu ), M Hd = m (l + 6 Hd ), 

NU3 : M q3 = mo(l + S q3 ), M u3,d3 = ™o(l + <W ( 2 - 2 ) 

NUG : Mi = mi/2, M 2 , 3 = m 1/2 (l + <5m 2)3 )- 

In the above 5h u and 5ff d define the nonuniversalities for the up and down Higgs mass 
parameters, M q3 is the left-handed squark mass for the 3rd generation, and M u3 (M^) are 
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the right-handed u-squark (d-squark) masses for the 3rd generation. The nonuniversalities in 
the gaugino sector are parameterized here by 5m 2 an d $m 3 - We have carried out a Monte 
Carlo scan with fiat priors using 10 6 model points in each of the three types of NUSUGRA 
models, taking the same input parameter ranges as specified in Eq,( p.l| ) and —0.9 ^ S ^ 1. 
Almost all of the mSP patterns seen for the mSUGRA cases were found in supergravity models 
with nonuniversal soft breaking, as the mSUGRA model is contained within the nonuniversal 
supergravity models. In addition we find many new patterns labeled NUSPs (nonuniversal 
SUGRA pattern), and they are exhibited in Table (||). As in the mSUGRA case one finds 
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Table 3: New 4 sparticle mass patterns that arise in NUSUGRA over and above the mSP patterns 
of Table (1). These are labeled nonuniversal SUGRA patterns (NUSP) and at least 15 new patterns 
are seen to emerge which are denoted by NUSP1-NUSP15. 

several pattern classes, CPs, SUPS, SOPs, and HPs as exhibited below. In addition, we find 
several Gluino patterns (GP) where the gluino is the NLSP. 

1. Chargino patterns (CP) : NUSP1, NUSP2, NUSP3, NUSP4 

2. Stau patterns (SUP) : NUSP5, NUSP6, NUSP7, NUSP8, NUSP9 

3. Stop patterns (SOP) : NUSP10, NUSP11 

4. Higgs patterns (HP) : NUSP12 

5. Gluino patterns (GP) : NUSP13, NUSP14, NUSP15. 

It is interesting to note that for the 4 sparticle landscape we find saturation in the number of 
mass hierarchies that are present. For example, for the case \i > in mSUGRA , increasing 
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the soft parameter scan from 1 x 10 parameter model points to 2 x 10 model points does not 
increase the number of 4 sparticle patterns. In this context it becomes relevant to examine as 
to what degree the relic density and other experimental constraints play a role in constraining 
the parameter space and thus reducing the number of patterns. This is exhibited in Table (|j) 
where we demonstrate how the relic density and the other experimental constraints decrease 
the number of admissible model points in the allowed parameter space for the mSUGRA 
models with both [i > and \i < 0, and also for the cases with nonuniversalities in the 
Higgs sector, nonuniversalities in the third generation sector, and with nonuniversalities in 
the gaugino sector. In each case we start with 10 6 model points at the GUT scale, and find 
that the electroweak symmetry breaking constraints reduce the number of viable models to 
about 1/4 of what we started with. We find that the allowed number of models translates into 
SUGRA mass patterns which are typically less than 100. The admissible set of parameter 
points reduces drastically when the relic density constraints are imposed and are then found 
to typically reduce the number of models by a factor of about 200 or more, with a reduction in 
the number of allowed patterns by a factor of 2 or more. Inclusion of all other experimental 
constraints further reduces the number of admissible points by a factor between 30% and 
50%, with a corresponding reduction in the number of patterns by up to 40%. The above 
analysis shows that there is an enormous reduction in the number of admissible models and the 
corresponding number of hierarchical mass patterns after the constraints of radiative breaking 
of the electroweak symmetry, relic density constraints, and other experimental constraints are 
imposed. 
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Trial 


Output 


No. of 


Relic Density 


No. of 


All 


No. of 


Type 


Models 


Models 


Patterns 


Constraints 


Patterns 


Constraints 


Patterns 


mSUGRAO > 0) 


10 6 


265,875 


55 


1,360 


22 


902 


16 


mSUGRAO < 0) 


10 B 


226,991 


63 


1,000 


31 


487 


18 


NUHO > 0) 


10 6 


222,023 


59 


1,024 


24 


724 


15 


NU3(m > 0) 


10 s 


229,928 


73 


970 


28 


650 


20 


NUGO > 0) 


10 6 


273,846 


103 


1,788 


36 


1,294 


28 



Table 4: An analysis of mass patterns for the four lightest sparticles. Exhibited in the table are 
the model type, the number of trial input points for each model, the number surviving the radiative 
electroweak symmetry breaking scheme as given by SuSpect (column 3), the number surviving when 
the relic density constraints are applied with MicrOMEGAs (column 5), the number surviving with 
inclusion of all experimental collider constraints (column 7), along with the corresponding number of 
hierarchical mass patterns in each case (column 8). 



2.3 Hierarchical patterns for the full sparticle spectrum 

We discuss now the number of hierarchical mass patterns for the full set of 32 sparticles 
in SUGRA models when the constraints of electroweak symmetry, relic density, and other 
experimental constraints are imposed. The result of the analysis is given in Fig.(|2|) and Table 
(H). Here one finds that increasing the number of model points in the scan does increase 
the number of patterns. However, the ratio of the number of patterns to the total number of 
models that survive all the constraints from the scan decreases sharply as shown in the right 
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2000 



Number of mSUGRA Models Number of mSUGRA Models 

Figure 2: Left panel: The number of hierarchical mass patterns for 32 sparticles vs the number of 
trial points for mSUGRA models which survive the electroweak symmetry breaking constraints, the 
relic density and all other experimental constraints. The number of hierarchical mass patterns show 
a trend towards saturation. Right panel : A similar phenomenon is seen in the ratio between the 
number of patterns over the number of surviving trial points in mSUGRA models. 



Models [No.] 


No. after constraints 


No. of patterns 


mSUGRA(/x > 0) [10 6 ] 


902 


505 


mSUGRA(^ < 0) [10 6 ] 


487 


268 


NUH(/i > 0) [10 6 ] 


724 


517 


NU3(^ > 0) [10 6 ] 


650 


528 


NUGO > 0) [10 6 ] 


1294 


1092 


All Above[5 x 10 B ] 


4057 


2557 



Table 5: The table exhibits a dramatic reduction of the landscape from upward of ~ O(10 28 ) hier- 
archical mass patterns for the 32 sparticle masses to a much smaller number when the electroweak 
symmetry breaking constraints, the relic density constraints, and other experimental constraints are 
applied. Column 1 shows one million input parameter points for each of the models investigated, and 
the number surviving all the constraints are exhibited in column 2, while column 3 gives the number 
of hierarchical patterns. 



panel of Fig.(^). This means that although saturation is not yet achieved one is moving fast 
towards achieving saturation with a relatively small number of allowed patterns for all the 32 
sparticles within SUGRA models consistent with the various experimental constraints. The 
analysis of Table (||) shows that the number of allowed patterns for the 32 sparticles, which 
in the MSSM without the SUGRA framework can be as large as O(10 28 ) or larger, reduces 
rather drastically when various constraints are applied in supergravity models. We note that 
some patterns are repeated as we move across different model types listed in the first column 
of Table (|5|). Thus the total number of patterns listed at the bottom of the last column of 
this table is smaller than the sum of patterns listed above in that column. We note that the 
precise number and nature of the patterns are dependent on the input parameters such as 
the top mass and a significant shift in the input values could modify the pattern structure. 



- 11 - 



60 



50 



40 



30 



20 



10 




ft 4 ft w ♦ ♦ ♦ 

♦ « ♦ * 'o^o **t ** » * * S ♦ 



-So 



10 



> 
u 

b 2 

o 

E 



5 O 0° . 

o o 8 oo 

- *° A 



<= „ o o o 6 ° 

° ° O + ° 



* „ ° ft ° % 

ft c 

o ° 0° ° ° 

4 o »o, °» ° 

o o o 

,+ 0° ° , 



J4« 

or™ 



o O 

o ° ♦o'co 

°* ^ ' 

1 o£o#° t Jtt 
o 0*00 J, <» O?" 




m„, (TeV) 



50 1 , 

o mSP1 

o mSP2 

o mSP3 

o mSP4 

40 ♦ mSP5 

♦ mSP6 

♦ mSP7 

♦ mSP8 

♦ mSP9 
30 4 mSP11 

4 mSP12 

4 mSP13 

» mSP17 

» mSP18 

20 » mSP19 

» mSP20 

» mSP21 

♦ mSP22 




mSUGRA 

\i<0 



4 . o O ! 

4 4 f Oq 



* « * 



10. « 



» * • 4% O " ' \,*' ''ft* . 



-10 



10 



3 



> 

to 

b 2 

o 

E 



o°8 ° 06 £>' 

O O o OO * ft 



% * 



O O fp _ "go 
o 

otfb 



jJ?JtlW** mSUGRA, n<0 



0.5 1 

m, /2 (TeV) 



1.5 



Figure 3: The dispersion of mSPs arising in mSUGRA in the tan/3 vs Aq/itiq plane (left panels), and 
in the mo vs m\/2 plane (right panels) for the /i > case (upper panels) and /i < case (lower panels). 
The analysis is based on a scan of 10 6 trial model points with flat priors in the ranges mo < 4 TeV, 
m \/2 < 2 TeV, 1 < tan/3 < 60, and |Ao/mo| < 10. mSPl is confined to the region where |Ao/mo| < 2. 
For the case \i < 0, no HPs are seen, and also, no model points survive in the region where tan/3 > 50 
in contrast to the jj. > case where there is a significant number for tan/3 > 45. 



3. Sparticle Patterns and the Nature of Soft Breaking 

3.1 Correlating mass hierarchies with the soft parameter space 

It is interesting to ask if the patterns can be traced back to some specific regions of the 
parameter of soft breaking from where they originate. This indeed is the case, at least, for 
some of the patterns. The analysis illustrating the origin of the patterns in the parameter 
space is given in Fig.(||). Exhibited are the landscape of sparticle mass spectra in the planes 
(I) tan/3 vs Aq/ttlq and (II) ttiq vs m^, when the soft parameters are allowed to vary in 
the ranges given in Eq.([2.1|). Many interesting observations can be made from these spectral 
decompositions. For example, a significant set of the mSPl (CP) models lie in the region 
lA)/ 77 ^)! < 2 and correspond to the Hyperbolic branch(HB)/Focus Point (FP) Q regions, 
while most of the SOPs have a rather large ratio of Aq/tjiq with the satisfaction of REWSB. 
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Figure 4: Dispersion of patterns in the tuq vs to 1/2 plane for fixed values of tan/3 and Ao/mQ. The 
region scanned is in the range mo < 4 TeV and < 2 TeV with a 10 GeV increment for each 

mass. Only a subset of the allowed parameter points relative to Fig.(^) remain, since the scans are on 
constrained surfaces in the mSUGRA parameter space. 

In this analysis we require that there be no charge or color breaking (CCB)|64], ^] at the 
electroweak scale. We note in passing that it has been argued that even if the true minimum 
is not color or charge preserving, the early universe is likely to occupy the CCB preserving 
minimum and such minima may still be acceptable if the tunneling lifetime from the false to 
the true vacuum is much greater than the present age of the universe^]. Next, we note that 
for the mSUGRA [i > case, the region around tan (3 = 50 has a large number of models 
that can be realized, while the region around tan (3 = 30 has far less model points. We also 
note that most of the HPs reside only in the very high tan/3 region in mSUGRA, but this 
situation can be changed significantly in the NUH case where HP points can be realized in 
the tan/3 region as low as tan/3 ~ 20. In the tuq vs mj/2 plane, one finds that most of CPs 
and HPs have a larger universal scalar mass than most of the SUPs and SOPs. 

Often in the literature one limits the analysis by fixing specific values of Aq and tan (3. 
For Aq the value most investigated is A$ = 0. However, constraining the values of Aq or of 
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tan (3 artificially eliminates a very significant part of the allowed parameter space where all the 
relevant constraints (the REWSB constraint as well as the relic density and the experimental 
constraints) can be satisfied as seen in Fig.(^). One can extract the familiar plots one finds 
in the literature where A$ and tan (5 are constrained from a reduction of the top-right panel 
of Fig.(^). The results of this reduction are shown in Fig.(Q) with a focused scan in specific 
regions of the soft parameter space. Specifically the bottom-left and top-right panels of Fig. (Q) 
show the familiar stau coannihilation |37], |68| , [5l[l regions and the HB/FP branch, the bottom- 
right panel gives the stau coannihilation region and the stop coannihilation region because 
of the relatively large Aq value, and the top-left panel is of the form seen in the works of 
Djouadi et al. |^3| where the Higgs funnel plays an important role in the satisfaction of the 
relic density. 
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Figure 5: An exhibition of the NUSPs and mSPs for the NUH, NU3, and NUG models in the tan/3 
vs Ao/mo plane. The range of SUGRA parameters are the same as the case mSUGRA (/i > 0). One 
may notice that the mSPl points arising from NU models lie in a relatively larger Ao/mo region. Most 
of the models in NU cases are still mSPs, and among the NUSPs, only two patterns have a relatively 
large population, these being NUSP1 and NUSP13. One may also notice that in NUH case, the HPs 
can exist in a low tan/3 region as opposed to the mSUGRA case where HPs can either exist in the 
large tan/3 region (/x > 0) or are totally eliminated (fj. < 0). 
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An analysis similar to that of Fig.(^) for the nonuniversal case is given in Fig. (|). Here 
in addition to the mSPs new patterns emerge which we label as nonuniversal sugra patterns 
or NUSPs. Among the NUSPs the dominant patterns are NUSP1 (CP) and NUSP13 (GP), 
which are seen to arise the model with nonuniversalities in the gaugino sector, i.e., the NUG 
model. In general, the NUG is dominated by the CP patterns whereas the NUH case is rather 
diverse offering the possibility of Higgs patterns at lower, less fine tuned values of tan (3. 

3.2 Benchmarks for sparticle patterns 

As discussed in Sec. (|2.lD , many of the sparticle mass patterns discussed in this analysis 
do not appear in the Snowmass, Post-WMAP, and CMS benchmark points. With some of 
these mSP and NUSP having a significant probability of occurrence, we therefore provide a 
larger set of benchmark points for the various patterns in different SUGRA scenarios. These 
benchmark points are exhibited in Tables (Pj9|,10,n|JT^) of the Appendix. Each of these 



benchmarks satisfies the relic density and other experimental constraints with SuSpect linked 
to MicrOMEGAs. We have explicitly checked that the first mSP benchmark point in each of 
the tables can be reproduced by using SPheno, and SOFTSUSY by allowing minor variations 
on the input parameters. The benchmarks are chosen to cover wide parts of the SUGRA 
parameter space. We give these benchmarks, several for each mass pattern, as the search 
for SUSY from the point of view of mass patterns has important consequences for LHC 
experimental searches. Some of the patterns are correlated with certain well investigated 
phenomena such as the HB /FP branches of REWSB and the stau-neutralino co-annihilation 
regions. However, many of the patterns arise from multiple annihilation processes. 

4. LHC Signatures for Mass Patterns 

4.1 Event generation and detector simulation 

Before moving to the discussion of the LHC signatures arising from various mSPs and NUSPs, 
we first give a detailed description of our LHC simulation procedure. 

After the imposition of all the constraints mentioned in the previous sections, such as 
the relic density constraints from WMAP data, the constraints on the FCNCs, as well as 
mass limits on the sparticle spectrum, we are left with the candidate model points for the 
signature analysis. For each of these model points, a SUSY Les Houches Accord (SLHA) file 



@ is interfaced to PYTHIA 6.4.11 (7| through PGS4 f| for the computation of SUSY 
production cross sections and branching fractions. In this analysis, for signals, we have 
generated all of PYTHIA's 2^2 SUSY production modes using MSEL = 39 2 . Leading order 



cross sections from PYTHIA and leading order cross sections from PROSPINO 2.0 [73] were 



cross checked against one another for consistency over several regions of the soft parameter 



space. TAUOLA 74] is called by PGS4 for the calculation of tau branching fractions as 



2 More specifically this choice generates 91 SUSY production modes including gaugino, squark, slepton, and 
SUSY Higgs pair production but leaves out singly produced Higgs production. For further details, see jn|. A 
treatment of singly produced Higgs production in the context of sparticle mass hierarchies was included in the 
analysis of Ref . M . 
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controlled in the PYTHIA parameter card (.pyt) file. With PGS4 we use the Level 1 (LI) 



triggers based on the Compact Muon Solenoid detector (CMS) specifications |75], 58] and the 
LHC detector card. Muon isolation is controlled by employing the cleaning script in PGS4. 
We take the experimental nomenclature of lepton being defined only as electron or muon and 
thus distinguish electrons and muons from tau leptons. SM backgrounds have been generated 
with QCD multi-jet production due to light quark flavors, heavy flavor jets (bb, tt), Drell- 
Yan, single Z/W production in association with quarks and gluons (Z+ jets / W+ jets), and 
ZZ, WZ, WW pair production resulting in multi-leptonic backgrounds. Extraction of final 
state particles from the PGS4 event record is accomplished with a code SMART ( = SUSY 
Matrix Routine) written by us which provides an optimized processing of PGS4 event data 
files. The standard criteria for the discovery limit of new signals is that the SUSY signals 
should exceed either 5\/Nsm or 10 whichever is larger, i.e., Nsusy > MaxjS^/AsM, 10} and 
such a criteria is imposed where relevant. We have also cross checked various results of our 



analysis with three CMS notes |F(], 77, [78[ and we have found agreement with these works 
using SMART and PGS4 for signal and backgrounds. 

We note that several works where sparticle signatures are discussed have appeared re- 
cently 7£, |8(], 32, S3, 84]. However, the issue of hierarchical mass patterns and the 
correlation of signatures with such patterns has not been discussed which is what the analysis 
of this work investigates. 

4.2 Post trigger level cuts and LHC signatures 

Generally speaking, there are two kinds of LHC signatures: (i) event counting signatures, and 
(ii) kinematical signatures. We have investigated both of these for the purpose of discrimi- 
nating the sparticle mass patterns. We list our event counting signatures in Table (|6|), where 
we have carried out analyses of a large set of lepton + jet signals. In our counting procedure, 
only electron and muon are counted as leptons, while tau jets are counted independently. For 
clarity, from here on, our use of 'jet(s)' will exclude tau jets. Thus, for jet identification, we 
divide jets into two categories: b-tagged jets and jets without b-tagging, which we simply 
label as b-jets and non-b-jets (see also (7j|). There are some counting signatures that only 
concern one class of measurable events, for example, the number of events containing one 
tagged b-jet and any other final state particles. There are also types of signatures of final 
state particles with combinations of two or three different species. For instance, one such 
example would be the number of events in which there is a single lepton and a single tau. 

When performing the analysis of event counting, for each SUGRA model point, we impose 
global post trigger cuts to analyze most of our PGS4 data. Below we give our default post 
trigger cuts which are used throughout the paper unless stated otherwise. 

1. In an event, we only select photons, electrons, and muons that have transverse momen- 
tum > 10 GeV and \rf\ < 2.4, p = (j,e,/j). 

2. Taus which satisfy Pf. > 10 GeV and \rf\ < 2.0 are selected. 

3. For hadronic jets, only those satisfying Pj, > 60 GeV and \r} 3 \ < 3 are selected. 
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4. We require a large amount of missing transverse momentum, P™ lss > 200 GeV. 

5. There are at least two jets that satisfy the Pt and r\ cuts. 

Our default post trigger level cuts are standard and are designed to suppress the Standard 
Model background, and highlight the SUSY events over a broad class of models. 
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Table 6: The tables give a list of 40 counting signatures along with the kinematical signatures analyzed 
for each point in the SUGRA model parameter space. L = e,fi signifies only electrons and muons. 



The different kinematical signatures we investigated for the purpose of discriminating 
among sparticle mass patterns are also exhibited in Table (||). One may further divide the 
kinematical signatures into two classes: namely those involving transverse momentum Pt and 
those which involve invariant mass. For those involving Pt, we have investigated missing Pt 
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distributions and the effective mass, the latter being the sum of missing Pj- and Pt of all jets 
contained within an event. For the kinematical variables using invariant mass, we reconstruct 
such quantities for four different cases, i.e., the invariant mass for all jets, for e + e~ pair, for 
pair, and for t + t - pair. The reconstruction of the invariant mass of t + t~ pair is based 
on hadronically decaying taus (for recent analyses see |51||). 

4.3 Discrimination among mSPs in mSUGRA 

We turn now to a discussion of how one may distinguish among different patterns. The 
analysis begins by considering the 902 model points that survive our mSUGRA scan with 
10 6 trial points, and simulating their LHC signals with PGS4 using, for illustration, 10 fb _1 
of integrated luminosity at the LHC. In our analysis we will focus mostly on the counting 
signatures. Here the most useful counting signature is the total number of SUSY events after 
trigger level cuts and post trigger level cuts are imposed. All other counting signatures are 
normalized with respect to the total number of SUSY events passing the cuts and thus appear 
as fractions lying between (0,1) in our figures. To keep the analysis statistically significant, 
we admit only those points in the parameter space that generate at least 500 total SUSY 
events. 

We give now the details of the analysis. In Fig.(||), we investigate the signature space 
spanned by a variety of signature channels. The top left panel gives a plot with one signature 
consisting of events with one lepton and the second signature consisting of events with no 
leptons. It is seen that the stop patterns (SOPs) that survive the cuts are confined in a small 
region at the right-bottom corner and have a significant separation from all other mSPs. The 
panel illustrates the negligible leptonic content in stop decays. The top-right panel is a plot 
between two signatures where one signature contains a tagged b-jet while the other signature 
has no tagged b-jets. In this case one finds a significant separation of the CPs and HPs from 
SUPs and SOPs. The lower-left panel gives a plot where one signature has two tagged b-jets 
and the other signature has only one tagged b-jet. One again finds that the CPs and HPs 
are well separated from the SOPs and the SUPs for much the same reason as in upper-right 
panel. Finally, a plot is given in the lower-right panel where one signature is the average 
missing Pj< while the other signature involves events with no tagged b-jets. Again in this plot 
the CPs (which include mSP4) and HPs are well separated from the SOPs and SUPs. 

The analysis of Fig.(|6|) exhibits that for some cases, e.g., for the patterns CP and HP in 
the upper right hand corner of Fig.(|6|), the separation between the SUGRA prediction and 
the Standard Model background is strikingly clear, allowing for the identification not only of 
new physics but also of the nature of the pattern that leads to such a signature. 

We discuss now the possibility of discriminating sub-patterns within a given pattern class. 
An analysis illustrating this possibility is given in Fig. (0). Here the top two panels illustrate 
how the sub-patterns mSPl, mSP2, mSP4 within the chargino class (CP) are distinguishable 
with appropriate choice of the signatures. A similar analysis regarding the discrimination for 
the sub-patterns in the stau class (SUP) is given in the two middle panels. The lower-left 
panel gives an analysis of how one may discriminate the stop sub-patterns mSPll, mSP12, 
mSP13 in the stop class (SOP), and finally the lower-right panel shows the plots that allows 



- 18 - 




0.2 



z 

2 



0.1 



o mSP1 

o mSP2 

o mSP3 

o mSP4 

♦ mSP5 

♦ mSP6 

♦ mSP7 

♦ mSP8 

♦ mSP11 
» mSP12 
» mSP13 
+ mSP14 
+ mSP16 

♦ SM . SUP 

1 0.2 



o° + .° 

oo+ + + 




HP 



0.4 



0.3 



0.2 



0.1 




0.3 



800 



- 600 



£ 400 



0.5 



0.7 



0.9 



Ob/N 



0.3 



0.4 



20( 



o mSP1 




♦ ♦ 


o mSP2 






> mSP3 




: 


o mSP4 






» mSP5 
♦ mSP6 






- ♦ mSP7 
♦ mSP8 
» mSP11 
mSP12 


+ 

+ 

HP *°° 




» mSP13 
+ mSP14 


++ + 

'o=V4 8 OO 

o ego* 


V ♦ a$ sop 


+ mSP16 




_ • SM 


* V 4 .\ 








mSP4 




o 






CP 


(^*) «SM 



1b/N 



0.3 



0.5 
Ob/N 



0.7 



0.9 



Figure 6: Top Left: An exhibition of the mSPs in the 1L vs 0L where the fraction of events to the total 
number of events in each case is plotted. The analysis shows that the Stop Patterns (SOP) appearing 
on the right-bottom corner are easily distinguished from other patterns. The analysis shows that SOP 
has few lepton signals. Top Right and Bottom Left: Plots in the signature space with fraction of 
events with lb vs 0b and 2b vs lb exhibiting the separation of CPs and HPs from SOPs and SUPs, 
with CPs and HPs occupying one region, and SOPs and SUPs occupy another in this signature space 
except for a very small overlap. Bottom Right: An exhibition of the mSPs in the signature space with 
the average missing for each parameter point in the mSUGRA parameter space along the y-axis 
and the fraction of events with 0b along the x-axis. The plot shows a separation of the CPs and HPs 
from SOPs and SUPs. Further, mSP4 appears isolated in this plot. Most of the CPs and HPs have 
less than 60% events without b-jet content. The ratios for the SUSY models refer to the SUSY signal 
only. The SM point is purely background. 



one to discriminate the Higgs patterns mSP14 and mSP16 from each other. There are a 
variety of other plots which allow one to discriminate among patterns. With 40 counting 
signatures one can have 780 such plots and it is not possible to display all of them. A global 
analysis where the signatures are simultaneously considered for a large collection of mSPs 
and NUSPs is discussed in Sec. (4.7). 

As mentioned in the above analysis we have included models which can produce at least 
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Figure 7: An exhibition of how the mSPs can be discriminated within a given class, i.e., within CPs, 
SUPs, SOPs, and HPs. The analysis shows that patterns within a given class can be discriminated. 



500 SUSY events with 10 fb which is lower than our estimated discovery limits for total 
SUSY events which are about 2200 in this case. The reason for inclusion of points below 
the discovery limit in the total SUSY events is that some of them can be detected in other 
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channels such as in the trileptonic channel while others will be detectable as the luminosity 
goes higher. We note in passing that reduction of admissible points makes separation of 
patterns easier. 



4.4 Sparticle signatures including nonuniversalities 
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Figure 8: Discrimination among mSPs within both mSUGRA and NUSUGRA models. Two mSPs 
are presented in each figure in different signature spaces to show the separation for each case. Signals 
are simulated with constant number of events in PGS4 for each pattern. 



In this subsection, we give an analysis including nonuniversalities in three different sectors: 
NUH, NU3, and NUG. In our analysis we simulate various models with the same constant 
number of events N which we take as an example to be N= I0 4 . To discriminate among 
the patterns in the signature space, we introduce another set of post trigger cuts, which we 
denote as 'b jet cuts', in addition to the default post trigger cuts specified in Sec. ( [4. 2D . The 
criteria in the b-jet cuts are the same as the default post trigger cuts, except that we change 
the condition 'at least two hadronic jets in the event' to 'specifically at least one b-tagged 
jet in the event'. We exhibit our analysis utilizing both the default cuts and the b jet cuts 
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in Fig.(||). One can see that even with inclusion of a variety of soft breaking scenarios, some 
mSPs still have very distinct signatures in some specific channels. 

Thus in the top-left panel we give a plot of mSP7 (SUP) and mSPll (SOP) in the 
signature space 1L/N (b jet cuts) vs OLlb/N, where OLlb/N is obtained with the default 
post trigger cuts. Here we find that these two model types are clearly distinguishable as 
highlighted by shaded and unshaded regions. A similar analysis with signatures consisting 
of lLlb/N (b jet cuts) vs OLlb/N for mSP4 (CP) and mSP7 (SUP) is given in the top- 
right panel. The lower-left panel gives an analysis of mSP4 (CP) and mSP5 (SUP) also in 
the signature space consisting of lLlb/N (b jet cuts) vs OLlb/N. Finally, in the lower-right 
panel we give an analysis of mSP3 (CP) and mSPll (SOP) in the signature plane e + /N vs 
lb/N. These analyses illustrate that the patterns and often even the sub-patterns can be 
discriminated with the appropriate choice of signatures for a general class of SUGRA models 
including nonuniversalities . 

4.5 The trileptonic signal as a pattern discriminant 
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Figure 9: A plot of the number of trilepton events versus the light chargino mass for three patterns, 
one from each class, CP, SUP and HP. The SUP pattern gives the largest trileptonic signal followed 
by the HP and CP patterns. 



The trileptonic signal is an important signal for the discovery of supersymmetry. For on- 
shell decays the trileptonic signal was discussed in the early days in [11, 85] and for off-shell 
decays in |$6| . (For a recent application see j78| ) . Here we discuss the trileptonic signal in the 
context of discrimination of hierarchical patterns. In Fig. @ we exhibit the dependency of the 
trilepton signal on the chargino mass. It is seen that mSP5 gives the largest number of events 
in this channel while the CP pattern (mSPl) and the HP pattern (mSP14) can also produce 
a large number of trilepton events above the discovery limit, while the chargino mass reach 
is extended for the mSP5 as opposed to the mSPl and mSP14. The above observations hold 
for some of the other SUP patterns as well. Thus the trileptonic signal is strong enough to 
be probed up to chargino masses of about 600 GeV in the SUP pattern. Another interesting 
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Figure 10: The number of tri-lepton events versus the sparticle mass splittings. The left panel 
shows clear separations for hierarchical mass patterns in the number of trilepton events produced with 
10 fb _1 as a function of the NLSP and the LSP mass splitting for the chargino (CP) pattern mSPl 
and Stau (SUP) mSP5. The plot on the right shows a similar effect for the case where the mass 
splitting is taken to be the difference of the CP odd Higgs boson mass and the LSP for both the Higgs 
pattern mSP14 and the stau pattern mSP5. The Standard Model background is highly suppressed in 
this channel. 



display of the trileptonic signal is when this signal is plotted against some relevant mass 
splittings. Thus the left-panel of Fig. (|i~0| ) gives an analysis for the trileptonic signal for two 
patterns: the Chargino pattern mSPl and the Stau pattern mSP5 plotted against the NLSP- 
LSP mass splitting with 10 fb _1 of data. The analysis of the left-panel of Fig. ( |l0|) shows 
that the SUP pattern presents an excellent opportunity for discovering SUSY through the 3 
lepton mode. The analysis also shows a clear separation among mass patterns and further a 
majority of the model points stand above the discovery limit which in this channel is ~ 15 



events under the post trigger level cuts discussed in Sec. (|4.2D . The right-panel of Fig. (|l0|) 
gives an analysis of the trileptonic signal vs the mass splitting of the CP odd Higgs and the 
lightest neutralino LSP for patterns mSP5 and mSP14. Again, we see a clear separation of 
model points. We note that CP odd Higgs can sometimes be even lighter than the LSP, and 
thus the quantity AM = Ma — M^.o plotted on the x-axis can sometimes become negative. 

4.6 Kinematical distributions 

In addition to the event counting signatures discussed above, the kinematical signatures are an 
important tool for pattern discrimination. We illustrate this using the kinematical variables 
consisting of missing Pt and the effective mass (see Table (||) for their definitions) and an 



illustration is given in Fig. (|Tl]) . Specifically the analysis of Fig. (11) uses four mSUGRA 
points one each in the patterns CP, SUP, SOP and HP . The analysis of Fig. (11) shows that 
the distributions for the CP, HP, SOP and SUP are substantially different. It is interesting to 
note that in the missing Pt distribution, the HP and SUP model points have a relatively flat 



- 23 - 




Figure 11: An exhibition of the missing Pt and of the effective mass distributions for 4 different 
mSUGRA models with each corresponding to one class of mSPs, and for the Standard Model. In 
the missing Pt distribution as well as in the effective mass distribution, the Standard Model tends 
to produce events with a lower missing Pt and a lower effective mass relative to the mSUGRA case 
which generates events at relatively higher missing Pt and effective mass. Further, there is a large 
variation between different mSUGRA models, as can be seen above. Thus, for example mSP5 (a stau 
pattern) and mSP14 (a Higgs pattern) have peaks at larger values of missing Pt and larger values of 
the effective mass relative to mSPl (a chargino pattern) and mSPll (a stop pattern). Additionally, 
the shapes of the distributions are also different. Only trigger level cuts are employed here. 



distribution compared to the CP and SOP model points. The missing Pt distribution and 
the effective mass distribution are useful when designing post trigger level cuts to optimize 
the signal over the background. For instance, one can take a 1 TeV effective mass cut to 
analyze the SUP and HP signals shown in Fig. (|TT|) , but this method will will not work well 
when it comes to the CP and SOP points since most of their events have a rather small 
effective mass. To illustrate that different models have different effective mass distributions, 
and consequently different effective mass cuts are needed for different patterns, an analysis is 
given in Fig.(^) for the same set of points in Fig.([iT|) with post trigger level cuts imposed. 

We also investigate the invariant mass distribution for the opposite sign same flavor 
(OSSF) di-leptons (e + e~, in Fig.([D|). We applied the default post trigger cuts as in 
Sec.([4.2|) to suppress the SM background. As a comparison the dominant Standard Model ti 



background is also exhibited. We have cross checked our work with the CMS Note [76], and 
found good agreement regarding the SUSY signals and the Standard Model background. It is 
seen that the two mSP points plotted in Fig. (|l3|) are easily distinguishable from each other. 

4.7 A 'global' analysis, fuzzy signature vectors, and pattern discrimination 

In the above we have given specific examples of how patterns can be differentiated from each 
other. In the previous sections we used only a few of the 40 signatures exhibited in Table (^). 
However, in the analysis we have carried out we have examined all of them. Thus for each 
parameter point we have analyzed 40 signatures. We now define correlations among these 



- 24 - 



500 



400 



c 

§ 300 
w 

o 

ai 

f 200 

3 
Z 

100 



°6 500 1000 1500 2000 2500 3000 
Effective Mass (GeV) 

Figure 12: The effective mass distributions for 4 different mSUGRA models with each corresponding 
to one class of mSPs, and for the Standard Model. Post trigger level cuts are imposed here. The bin 
size used here is 25 GeV. We exhibit the mSUGRA points used here in the order (mo, m\/2, Aq, tan (3, 
sign//): CP (3206.9, 285.3, -1319.8, 9.7, +1), SUP (92.6, 462.1, 352.2, 4.5, +1), SOP (2296.9, 625.0, 
-5254.9, 13.6, +1), and HP (756.8, 387.0, 1144.9, 56.5, +1). 
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Figure 13: A plot of the opposite sign same flavor (OSSF) di-lepton invariant mass distribution at 
LHC with 10 fb" 1 with the default post trigger cuts imposed for two different mSP points on top of 
the SM it background. The two mSPs, mSP4 (1674.9, 137.6, 1986.5, 18.6, +1) and mSP5 (84.4, 429.3, 
-263, 3.4, +1), are clearly distinguishable from each other in the distribution. The mSP4 model point 
shown here has recently been investigated |?7j in the context of helicity amplitudes as a discovery 
mechanism for supersymmetry. 

signatures. Thus consider an ordered set where the signatures are labeled Si, S2, ••, S40 and 
let the number of events in each signature be iVi, A^, .., A/40. Define a signature vector for a 
given point x a (a = 1, 2, ..,p) in the parameter space 




^ttbar background 
mSP4 (CP) 
> mSP5 (SUP) 



Post Trigger Level Cuts Imposed 




(4.1) 



where = Nf/N and N is the total number of SUSY events. As the parameter point x a 
varies over the allowed range within a given pattern it generates a signature vector where the 
elements trace out a given range. Thus for a pattern X one generates a fuzzy pattern vector 
A£ x so that 

A£ x = (A£f,A£ 2 x ,..,A£f ), (4.2) 

where is the range traced out by the element ^ as the parameter point x a moves in 
the allowed parameter space of the pattern X. What makes the vector A£ fuzzy is that its 
elements are not single numbers but a set which cover a range. We define now the inner 
product of two such fuzzy pattern vectors so that 

C XY ^(A£ x |A£ y ) = 0(l) (4.3) 

where the inner product is if the element A£^ and A£jf overlap for all i (i = 1, .,40), and 
1 if at least one of the elements of pattern X, A£^ does not overlap with A£^, the element 
for pattern Y. Therefore, if for two patterns X and Y one finds there is no overlap at least for 
one signature component A£j, then these two patterns can be distinguished in this specific 
signature and one obtains Cxy = 1- Otherwise Cxy = which means that all components 
of A£ and A£ Y have an overlap and cannot be distinguished under this critera. We can 
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Figure 14: A table exhibiting the discrimination of patterns using the criterion of Eq.( |4.3| ) where 
various signatures with both the default post trigger cuts and b jet cuts are utilized. If the element of 
i th row and j th column is 1, i.e., Cy = 1, one can distinguish the i th mass pattern from the j th one. 



generalize the above procedure for the signatures 

Cm = #> (*,j = l,...,40). (4.4) 
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Repeating the previous analysis, one can construct another fuzzy signature vector for pattern 
X as 

AC X = (ACf 2 ,..,A£$,..,AC& j40 ) (4.5) 

where the elements have a range corresponding to the range spanned by the soft parameters 
x a as they move over the parameter space specific to the pattern. Further, the definition of 
the inner product Eq. ([4.3D still holds for this new fuzzy signature vector. We have carried 
out a full signature analysis of such comparisons, using 40 different signatures, and their 



combinations as defined in Eq. (4.4) and Eq. ([4.5]). An illustration of the global analysis is 
given in Fig. (|H]) . The analysis shows that it is possible to often distinguish patterns using 
the criterion of Eq.(|4.3j). We note that the analyses exhibited in Fig.(||) are the special cases 
of the results in Fig.(|l4}). For instance, the clear separation between mSP7 and mSPll in 
the signature space shown in the top-left panel of Fig. (^j) gives the elements C45 = C54 = 1 
of Fig. (|i~4|). As emphasized already the analysis of Fig. fli~4|) is for illustrative purposes as we 
used a random sample of 22 patterns out of 37. Inclusion of each additional mass pattern 
brings in a significant set of model points which need to be simulated, and here one is limited 
by computing power. The full analysis including all the patterns can be implemented along 
similar lines with the necessary computing power. Finally we note that the analysis in Fig.(|i~4j) 
is done without statistical uncertainties. Inclusion of uncertainties in pattern analysis would 
certainly be worthwhile in a future work. 

5. Signature Degeneracies and Resolution of Soft Parameters 
5.1 Lifting signature degeneracies 

It may happen that two distinct points in the soft parameter space may lead to the same 
set of signatures for a given integrated luminosity within some predefined notion of indistin- 
guishability. Thus consider two parameter points A and B and define the 'pulls' in each of 
their signatures by 



Pi 



1 A Ri 

nf — nf\ 



^AB 



yJiSnfY + (5nfY + (5nf M ) 2 . (5.1) 



Here 5nf ~ \J n t ls the uncertainty in the signature events nf-, and we estimate the SM 
uncertainty as 5nf M ~ ^Jy{bnf + 5nf)/2. Here the parameter y parameterizes the effect of 
the SM events, and for the analysis in this section, we take y = 1. In other words, if the 
pulls in each of the signatures is less than 5, then the two SUGRA parameter space points 
are essentially indistinguishable in the signature space. In such a situation one could still 
distinguish model points either by including more signatures, or by an increase in luminosity. 
Thus, for example, inclusion of the Higgs production cross sections, B s — > u + u~ constraints, 
as well as the inclusion of neutralino proton scattering cross sections constraints tend to 
discriminate among the model parameter points as shown in Ref. 0. Here we point out 
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37 
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7003 
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38 


2b 


34 


37 


0.3 


29 


39 


1.0 
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1.4 
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1.4 


39 


3b 


1 


1 


0.0 


6 


2 


1.1 
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1.7 
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0.4 


40 


4b 


1 





0.7 





2 


1.1 


27 


20 


0.8 


19 


29 
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Table 7: An exhibition of lifting the degeneracy of two points in the mSUGRA parameter space using 
luminosity. Two pairs of points (A, B) and (A', B') arc indistinguishable under the 2 sigma criteria 
at 10 fb _1 luminosity (column 3-8), but can be clearly separated when the luminosity increases to 500 
fb _1 (column 9-14). The Standard Model uncertainty is estimated as Snf M — (Snf + 5nf)/2. 



that in some cases increasing the luminosity can allow one to lift the degeneracies enhancing 
a subset of signatures in one case relative to the other. For illustration we consider the 
following two sets of points in the pattern mSP5 in the mSUGRA parameter space in the 
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following order (mo, Wi/2> Aq, tan/3, sign//). 



Point A (192.6, 771.3, 1791.1, 8.8, +1), 

Point B (163.0, 761.3, -775.8, 4.7, +1); ^ ' ' 

Point A' (159.3, 732.3, -783.1, 5.6, +1), 

Point B' (163.5, 753.3, -918.2, 3.3, +1). ^ ' ' 

In Table (D we compare the pulls for the pairs of points (A, B) and (A', B') at an integrated 
luminosity of 10 fb -1 and 500 fb -1 . For points A and B, one finds that the pulls are all 
less than 2 for an integrated luminosity of 10 fb -1 . However, for an integrated luminosity 
of 500 fb" 1 , the pulls for signatures (6,7,8, 11,21,34) increase significantly and the pull for 
signature number 7 is in excess of 4.5 allowing one to discriminate between the two parameter 
points A and B. A very similar analysis is carried out for parameter points A' and B' . Here 
one finds that the signature (0, 1, 2, 3, 6, 12, 14, 32, 33, 36) receive a big boost as we go from 10 
fb -1 to 500 fb -1 , and the signatures (0, 1,2,6,36) give pulls greater than 5, with the largest 
pulls being in excess of 12, allowing one to discriminate between the parameter points A' and 
B' . We note the analysis ignores systematic errors and also does not consider an ensemble 
of simulations. Nonetheless it does illustrate the effects of moving from a low to a high 
LHC luminosity allowing one to discriminate some model pairs, which appear degenerate in 
the signature space at one luminosity, but can become distinct from each other at a larger 
luminosity. 

5.2 Resolving soft parameters using LHC data 

We discuss now the issue of how well we can resolve the points in the parameter space x a (a = 
l,..,p) for a given luminosity. Consider Eq.( |5.lD and set 5N = y/N, and parameterize the 
standard model uncertainty by SN SM = ^/y5N. Next we set the criterion for the resolution of 
two adjacent points in the SUGRA parameter space separated by Ax a so that the separation 
in the signature space satisfies 

- AN =5. (5.4) 
V2iV + yN K ' 

Since TV = cr susy (x a )Ci J iic, where <T susy is the cross section for the production of sparticles, and 
jClhc is the LHC integrated luminosity, the resolution achievable in the vicinity of SUGRA 
parameter point x a at that luminosity is given by 

^fe+d^ (5 . 5) 

In Fig. (|T5|) we give an illustration of the above when ttiq varies between 500 GeV and 2000 



GeV while my 2 = 500 GeV, A = 0,tan/3 = 30, and \i > 0. From Fig.(|l5|) one finds that the 
resolution in m-o strongly depends on the point in the parameter space and on the luminosity. 
Quite interestingly a resolution as small as a few GeV can be achieved for mo in the range 
500-1000 GeV with 1000 fb -1 of integrated luminosity. A similar analysis varying m-i/2 in 
the range 500-900 GeV for the case when mo = 500 GeV, Aq = 0,tan/3 = 30 and /x > 0, 
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Figure 15: An analysis showing the resolutions in too and m 1 / 2 that can be reached with 1000 fb^ 1 
of integrated luminosity under the REWSB constraints. The two left panels give the number of SUSY 
events vs too (top left panel) and vs TO1/2 (lower left panel) for 1000 fb -1 of integrated luminosity. 
The right panels give the resolutions in too (top right panel) and in mi/2 (lower right panel) using the 
left panels. 



shows that a resolution in m.1/2 as low as 1 GeV can be achieved with 1000 fb of integrated 
luminosity. 

6. Conclusion 

The minimal supersymmetric Standard Model has 32 sparticle masses. Since the soft breaking 
sector MSSM is arbitrary, one is led to a landscape of as many as 10 28 or more possibilities 
for the sparticle mass hierarchies. The number of possibilities is drastically reduced in well 
motivated models such as supergravity models, and one expects similar reductions to occur 
also in gauge and anomaly mediated models, and in string and brane models. In this work 
we have analyzed the mass hierarchies for the first four lightest sparticle (aside from the 
lightest Higgs boson) for supergravity models. Specifically, in Sec.(|2|) we analyzed the mass 
hierarchies for the mSUGRA model and for supergravity models with nonuniversalities in 
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the soft breaking in the Higgs sector, nonuniversalities in the soft breaking in the third 
generation sector, and nonuniversalities in the soft breaking in the gaugino sector. It is found 
that in each case only a small number of mass hierarchies or patterns survive the rigorous 
constraints of radiative breaking of the electroweak symmetry, relic density constraints on 
cold dark matter from the WMAP data, and other experimental constraints from colliders. 
These mass hierarchies can be conveniently put into different classes labeled by the sparticle 
which is next heavier after the LSP. For the SUGRA models we find six different classes: 
chargino patterns, stau patterns, stop patterns, Higgs patterns, neutralino patterns, and 
gluino patterns. Benchmarks for each of these patterns were given in Sec.@. In Sec.(Q) 
we discussed the techniques for the analysis of the signatures and the technical details on 
simulations of sparticle events. In this section we also discuss the backgrounds to the SUSY 
phenomena arising from the Standard Model processes. Additionally we discussed here the 
identification of patterns based on 40 event identification criteria listed in Fig. (14). It is found 
that these criteria allow one to discriminate among most of the patterns. An analysis of how 
one may lift degeneracies in the signature space, and how accurately one can determine the 
soft parameters using the LHC luminosities is given in Sec.(^). It is hoped that the analysis 
of the type discussed here would help not only in the search for supersymmetry but also allow 
one to use the signatures to extrapolate back to the underlying supersymmetric model using 
the experimental data when such data from the LHC comes in. In the above our analysis was 
focused on supergravity unified models. However, the techniques discussed here have a much 
wider applicability to other models, including models based on gauge and anomaly mediated 
breaking, as well as string and brane based models. 
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8. Appendix: Benchmarks for sparticle mass hierarchies 



Chargino Patterns (CPs) 
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Table 8: Benchmarks for the class CP where the chargino xt is the NLSP in mSUGRA and in 
NUSUGRA models. Benchmarks are computed with m b m {m b ) = 4.23 GeV, a s M (M z ) = .1172, and 
m t (pole) = 170.9 GeV with SuSpect 2.34 interfaced to micrOMEGAS 2.07. 



Gluino Patterns (GPs) 
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Table 9: Benchmarks for the class GP where the gluino g is the NLSP. Such a pattern was only seen 
to appear in NUSUGRA models with non universal gaugino masses. 
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Stau Patterns (SUPs) 
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Table 10: Benchmarks for the class SUP where the stau t\ is the NLSP in mSUGRA and in 
NUSUGRA. 
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Stop Patterns (SOPs) 
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Table 11: Benchmarks for the class SOP where the stop t\ is the NLSP in mSUGRA and in 
NUSUGRA models. 



Higgs Patterns (HPs) 
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Table 12: Benchmarks for the class HP where the Higgs boson (A, H) is the next nearest heavy 
particle after the LSP in mSUGRA and in NUSUGRA. The LSP and (A, H) sometimes are seen to 
switch. 
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